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Abstract. A notion of geometric symmetry is introduced that generalizes the 
classical concepts of Killing fields and other affine coUineations. There is a sense in 
which flows under these new vector fields minimize deformations of the connection 
near a specified observer. Any exact affine coUineations that may exist are special 
cases. The remaining vector fields can all be interpreted as analogs of Poincare and 
other well-known symmetries near timelike worldlines. Approximate conservation 
laws generated by these objects are discussed for both geodesies and extended 
matter distributions. One example is a generalized Komar integral that may be 
taken to define the linear and angular momenta of a spacetime volume as seen by 
a particular observer. This is evaluated explicitly for a gravitational plane wave 
spacetime. 



PACS numbers: 02.40.-k, 04.20.Cv, 04.25.-g 

1. Introduction 

It is rarely possible to model realistic physical systems with exact solutions to the 
equations of some general underlying theory. Despite this, many interesting problems 
deviate only slightly from a model problem that can be understood exactly. Such 
solutions are usually tractable only because of symmetry assumptions. Once they're 
understood, perturbation theory may be used to understand many systems that 
"almost" satisfy the appropriate symmetry principle. While this statement has a clear 
intuitive meaning, quantifying it can be difficult. It is also not necessarily obvious how 
- or if it is meaningful - to uniquely propagate a symmetry from solutions where it 
is exact into the perturbations that break it. These issues are particularly important 
in the context of conservation laws. As an example, one might want to know how to 
construct approximately conserved quantities that are unique generalizations of some 
exact counterpart in a similar system. There would hopefully be a sense in which such 
quantities varied slowly for some class of small perturbations. 

Some steps towards understanding problems like these are explored here in the 
context of affine coUineations (of which Killing vectors are special cases) associated 
with curved spacetimes. While these kinds of symmetries rarely exist, there are various 
senses in which approximate replacements can usually be introduced. One method for 
finding vector fields that are "almost Killing" is to write down an action whose value 
provides some sense for how nearly a particular flow preserves the metric [1^ . There are 
important caveats to this interpretation, although the final result is that any vector 
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field extremizing such an action satisfies a fairly simple generalization of Killing's 
equation. Various reasons have been given for suggesting other extensions as well 
[21 [21 m m H] ■ While any genuine Killing vectors that might exist are solutions to all of 
these equations, it is not usually clear how the remaining fields should be interpreted. 
This problem arises even in flat spacetime. 

What form an approximate symmetry should take is highly dependent on its 
intended use. One application is in the estimation of a black hole's angular momentum. 
This requires finding rotational Killing fields on certain 2-spheres foliating a horizon. 
Various methods have therefore been developed for defining such objects using only 
the intrinsic geometry of these surfaces [Zl [3 [9] . The concept of approximate Killing 
fields has also been adapted for use on initial data sets used in 3 + 1 splits of Einstein's 
equation [10]. 

The approach taken here is to define a set of vector fields in a four dimensional 
volume that can all be viewed as analogs of known symmetries in Minkowski spacetime. 
The physical interpretation is that these fields may be viewed as generators of 
approximate symmetries by a specified observer. Any sufficiently small region near 
a particular point can be made to look nearly flat. Some structures from Minkowski 
spacetime may be therefore be introduced very near this point. Approximate 
symmetries that take advantage of this fact are proposed in Sect. [2j It is then 
shown in Sect. [3| that analogous objects can also be introduced near an observer's 
worldline. These sorts of vector flelds can actually be extended in a non-perturbative 
way to finite regions around the point or worldline from which they were constructed. 
A well-defined subset provides a precise analog of the Poincare group. Translations, 
rotations, and boosts very near an observer extend in a useful way to cover large 
portions of the spacetime. Any exact symmetries that may exist are included as 
special cases. Some connections to conservation laws are discussed in Sect. [H and a 
simple example involving gravitational plane wave is finally presented in Sect. [5] 



Exact symmetries 

Given some spacetime {A4,gab), there are several types of exact symmetries that 
may be discussed. The most common of these take the form of vector fields whose 
associated diffeomorphisms preserve some geometric structure. The most ubiquitous 
examples are the Killing fields. Their flows preserve the metric. Vector fields with 
this property satisfy 

i^YKgab = 0. (1) 

Any solutions that may exist can be used to flnd conserved quantities associated with 
geodesies or matter distributions [TT], identify mass centers [T^], simplify Einstein's 
equation [13l [M] , classify its solutions [Ml [15] , and so on. 

This utility has (among other reasons) motivated various generalizations of ([T]). 
Perhaps the simplest of these arises from considering flows that preserve the metric 
only up to some constant factor: 

C-Yy^gab = 2cgQb. (2) 

Any satisfying this equation with constant c is known as a homothetic vector field 
or homothety. Allowing the dilation factor c to vary would define a conformal Killing 
vector. These objects preserve the metric up to an arbitrary multiplicative factor. 
The standard Killing vectors are special cases of either class. Like them, conformal 
and homothetic vector fields usually do not exist. Their presence can be very useful. 
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however. The existence of a proper (non-Killing) homothetic vector is often used to 
define a notion of geometric self-similarity, for example. Such objects therefore appear 
in certain models of gravitational collapse and cosmology. They are also related to 
the appearance of critical phenomena in general relativity |16j . 

A simple generalization of the homotheties can be found by considering vector 
fields that satisfy 

^a^YA9bc = 0. (3) 

Solutions to this equation are known as affine coUineations. They are the generators 
of infinitesimal affine transformations. Killing and homothetic vector fields are special 
cases. All of the affine coUineations may be interpreted geometrically as preserving 
the Levi-Civita connection. This means that Cy^ ^-nd commute when acting on 
arbitrary tensor fields. Geodesies and their affine parameters are also preserved under 
the action of any Y^. Although this might seem to be a significant generalization of the 
Killing vector concept, solutions rarely exist. The only non-flat vacuum spacetimes 
that admit non- homothetic affine coUineations are the pp-waves |17) . Similarly, it 
has been shown that proper homotheties cannot exist in any asymptotically flat 
vacuum spacetime with positive Bondi mass |18] . Despite these results, interesting 
affine coUineations can occasionally be identified in geometries that are not Ricci-flat. 
Doing so provides a number of simplifications for various problems. Some of these 
derive from the fact that 

Kab = Cvj^gab (4) 
is a second-rank Killing tensor; i.e. 

^(aKbc) = 0. (5) 

It should be noted that not all symmetric tensors satisfying this equation can be 
derived from affine coUineations. One counterexample is the Killing tensor associated 
with Carter constants in Kerr. 

Transformations generated by affine coUineations can be viewed as mapping 
geodesies into geodesies. They preserve the affine parameters of each curve. Dropping 
this latter requirement recovers the so-called projective coUineations. A precise 
definition may be found in [15], although it will not be needed here. One of their 
interesting consequences is that they leave invariant the projective curvature tensor: 

CY^{R\cd-\5l^Rd]b)=0. (6) 

Vector fields satisfying this equation are not always projective, however. 

The list of definitions here could keep growing as new fields are added that 
preserve more and more geometric structures. Interestingly, the quantities introduced 
so far all share a very useful characteristic that does not easily generalize: the space of 
vector fields in each of the mentioned classes has finite dimension. Furthermore, any 
single element is uniquely determined by its value and the values of its first one or two 
derivatives at a single point. These properties are well-known for Killing fields. Four 
dimensional spacetimes (which is all that will be considered here) admit a maximum 
of 10 linearly independent Killing vectors. At most one homothety can exist that 
is not itself Killing. The maximum number of (not necessarily proper) conformal 
Killing vectors is 15, and the affine coUineations total no more than 20. Finally, the 
space of projective coUineations has a maximum of 24 dimensions [15j . Properties like 
these do not hold for vector fields whose flows leave invariant the Riemann, Ricci, or 
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Einstein curvature tensors of a given spacetime. Despite this, the class of approximate 
symmetries introduced below is constructed so as to have finite dimension. Any given 
member is fixed by its value together with the values of its first derivatives at a point. 
Unlike the exact symmetries, these objects always exist at least in some finite region. 
After fixing a reference frame, the space of approximate symmetries has exactly 20 
dimensions. Ten of these will be identifiable as generalized Killing vectors, while the 
remaining ten will be related to more general affine coUineations. 

It has already been remarked that the presence of Killing fields implies the 
existence of various conserved quantities. The same can also be said for more general 
coUineations. Extensive discussions of exact symmetries and associated integrals of the 
geodesic equation may be found in [19] . Many of these symmetries are non-Noetherian 
in the sense that they preserve the equations of motion, but not the action. Despite 
this, their presence allows constants of motion to also be assigned to arbitrary stress- 
energy distributions satisfying Einstein's equation [501 [21] • As will be discussed in 
Sect, m generalizations of these quantities can be associated with any approximate 
affine coUineations that are identified. 

2. Symmetries near a point 

Generic symmetries in general relativity are usually discussed in the context of 
asymptotically flat spacetimes. There then exist approximate notions of isometry that 
improve as one approaches infinity Generalizations of these ideas also exist for 
geometries with somewhat more complicated (but still highly symmetric) asymptotic 
behavior like that of anti-de Sitter [52]. The assumption of a simple limiting form 
for the geometry makes it convenient to invariantly describe certain properties of a 
spacetime in terms of "measurements at infinity." Quantities that may be identified 
as a spacetime's total energy or angular momentum appear naturally, for example. 

While useful in many contexts, these ideas do not always translate into 
observations made by physical observers. Measurements like those expected from 
gravitational wave detectors do come very close to fitting into this formalism. Others 
can require a more local description. In particular, it is sometimes important to 
understand what given observers would experience inside strongly curved regions of 
spacetime. Abstracting the concept of an observer to a timelike worldline F, vector 
fields may be introduced in (say) some convex neighborhood ly of F that act like 
approximations to Killing fields or more general coUineations. This is always possible, 
and the symmetries these vectors generalize become exact on F itself. Limiting 
coUineations can evidently be useful on scales that are either very large or very small. 
It is much less clear how to easily describe systems at intermediate distances. 

2.1. Motivation 

The idea of a local symmetry just outlined is best introduced by first considering vector 
fields 'ip°'{x,"f) that generalize the affine coUineations in some reasonable way near a 
fixed reference point 7. Let these objects be defined inside a normal neighborhood N 
of this point. It is intuitively obvious that vector fields may always be chosen such 
that V aC^gbc vanishes at 7. While this condition is reasonable to require, it is not 
very interesting by itself. Much more can be said if each vector field in this class is 
uniquely fixed in N by knowledge of 

^'(7,7), Va^^(7,7)- (7) 
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It will be assumed that each 1/)° depends linearly on this initial data with no 
degeneracy. This implies that there are always 4 + 16 = 20 linearly independent 
vector fields defined about any given point in a four dimensional spacetime. Note 
that indices in ([7]) have been written in a sans-serif font to emphasize that they are 
associated with the preferred point 7. 

Approximate affine coUineations with the appropriate properties may be 
constructed by projecting symmetries of the tangent space T-^N into N using the 
exponential map. Consider the linear transformations 

X" -> X" + eBb"X^ (8) 

of vectors in this space parameterized by an arbitrary tensor B^^. Being a vector 
space, TjN has a preferred origin. Adding another term to ([5]) to shift that origin 
would be awkward. The translational symmetries that such a procedure might produce 
are certainly important, although generating them requires a more subtle treatment 
described below. For now, consider only the vector fields 

= X^Bb^ (9) 

associated with homogeneous transformations of the given form. These clearly satisfy 

^/:*..c = 2^(,,,(,)^v,^)=o, (10) 

so they are affine within T^N in the sense of ([3]). Such transformations can be made 
to induce shifts x x + eip in spacetime points associated with vectors X^ via 

x^exp^X. (11) 

A simple relation between ^/;° and is found by introducing Synge's world 
function a{x,y) — a{y,x). This two-point scalar returns one-half of the squared 
geodesic distance between its arguments. The assumption that iV be a normal 
neighborhood of 7 ensures that a(x, 7) is uniquely defined for all points x in this region. 
Many of its properties are reviewed in [23l [25 • Most importantly for the problem at 
hand, the first derivative of the world function effectively inverts the exponential map. 
Any set {7, x, X"-} satisfying (fTTjl is related via 

X,=-a,{x,j), (12) 

where the common shorthand da — Vgcr — da/d^^ has been used. The right-hand 
side of (fT^ generalizes the concept of a separation vector between two points. It is 
useful in that a straightforward expansion shows that linear transformations of the 
form ([5]) effectively shift spacetime points by an amount parameterized with a vector 

■!/;''(x,7) satisfying 

= -tr^V". (13) 

If the various components oi X^ as defined in are used as coordinates, the bitensor 
—a^a — —g^^d^cr/dx'^dj^ reduces to the identity. Components of tp"" and are 
therefore identical in normal coordinate systems of this type. In general, it is useful 
to introduce 

H% = [-<J\]-^ (14) 

as the matrix inverse of the operator appearing in (|13p . This always exists in the 
regions considered here. Using ^ now shows that 

^''(x,7) = -i^^abB^^ (15) 
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Holding 7 fixed, this equation defines a 16-paranieter family of vector fields generated 
by -Bab = VaV''b(7, 7)- Every such 7/)" vanishes at 7. It also satisfies ^ at this point. In 
flat spacetime, these vector fields coincide everywhere with exact affine collineations. 

Not all such symmetries are included in psp . however. The four translational 
Killing fields are missing. These can be obtained by considering transformations that 
directly shift the base point 7. Perturbations of this form cannot leave fixed, as 
the initial and final vectors must be elements of different spaces. Introducing some , 
we therefore demand that be parallel-transported along the curve that 7 follows 
under the one-parameter family of transformations 

j^j + eA. (16) 

Using this together with p2|) and the homogeneous transformation ([8]) generates the 
full 20-parameter family of approximate affine collineations 

^'^ ^H\{<j\A''-ai,B'''). (17) 

Given any A^ and B^^, these objects all satisfy 

Va/:^gbc(7) = 0. (18) 

The initial data 

A'=r{l,l), S^' = V^V^(7,7) (19) 
determine ip°'{x,j) throughout N. In Minkowski spacetime, one finds that 

= A" + (x - -ffBp" (20) 

in the usual coordinates. These coincide exactly with all of the affine collineations in 
this geometry. 

In general, vector fields satisfying ([3]) in a curved spacetime also have the form 
P7p for some A^ and B^^. This is most easily seen by noting that vector fields with the 
given form have been obtained before as general solutions to the equation of geodesic 
deviation (also known as the Jacobi equation) [25] 

aV^(VbV,^a - Rab/^d) = 0. (21) 
For any fixed x, this is an ordinary differential equation along the geodesic connecting 
that point to 7. Solving it repeatedly for all geodesies in N passing through this origin 
reproduces the vector fields P7|) . It is clear that such solutions always exist as long 
as the geometry is reasonably smooth. These are the spacetime's Jacobi fields about 
7. The bitensors iJ^aC^b and -ff^a^b are known as Jacobi propagators. 

Solutions to the geodesic deviation equation effectively map one geodesic into 
another while preserving the affine parameters of both curves. It was noted above that 
this is the defining characteristic of affine collineations. The difference is that such 
vector fields must map every geodesic into another geodesic. This intuitive argument 
makes it clear that all affine collineations ~ or Killing fields as special cases - must be 
solutions of (|2ip . The proof follows from noting that second derivatives of any exact 
affine collineation Y^{x) must satisfy 

VfcV.FX = ■'Rbdc'^Yl (22) 
This result is clear from (|A.12[) . and is actually equivalent to Substituting it 

into ((2T|) shows that all affine collineations are indeed special cases of Jacobi fields. 
As expected, satisfies the geodesic deviation equation along all geodesies; even 
those that do not pass through 7. This point illustrates precisely how (PTjl generalizes 
the equation defining an affine collineation. It is simply (|22p contracted into 
Alternatively, the Jacobi equation is equivalent to (jA.Sp . 
To summarize, the following is now evident: 
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Theorem 1 Let N be a normal neighborhood of some point 7. Define a Jacobi field 
4^°'{x,j) to be a solution of throughout this region. It is explicitly given by |J7[ ) 
for some initial data with the form il9\) . The set of all Jacobi fields about a fixed 7 
forms a 20- dimensional group in four spacetime dimensions. Each element satisfies 
-C^V = 7. Furthermore, all affine collineations are members of this group. 

These properties motivate our identification of the Jacobi fields as generalizations of 
affine collineations near 7. 

Further results that strengthen this decision are derived in the appendix. Even 
though few Jacobi fields are genuine afhne collineations, all can be interpreted as exact 
symmetries of certain quantities connected with the spacetime's geometric structure: 

Theorem 2 Given a Jacobi field defined as in theorem [II it is always true that 

Ctpcr^ ~ 'C^cr" = C^,a^a = ^4>H"^3 = 0, (23) 

where one of the arguments in each of these equations is taken to be the origin 7. 

Lie derivatives on two-point tensor fields are defined to act independently on each 
argument. See (jA.ip . for example. Quantities appearing in this theorem are all 
important in Riemann normal coordinate systems parameterizing arbitrary points 
X by the components of — —a^(x,"f). In terms of a more direct interpretation of 
the Jacobi fields as approximately satisfying ([3]), Lie derivatives of the metric with 
respect to an arbitrary ip'^ are strongly constrained by the identities (|A.7[) - (jA.10[) . 

Statements of this sort do not exhaust the connections between the Jacobi 
equation and a spacetime's symmetries. There is a sense in which higher-rank Killing 
tensors that may exist are also solutions to the geodesic deviation equation [2E\. 
Furthermore, projective collineations can be shown solve an inhomogeneous form of 
(PT|) proportional to cTa [57]. These observations will not be discussed any further here, 
although it is possible that they could be used to generalize the present framework. 

2.2. Special cases 

It is often useful to single out a subset of the Jacobi fields distinguished by 
antisymmetric Bat, — Va?Ab- These may be said to generalize only the Killing fields 
of a given spacetime. Distinguishing them with a subscript "K," they clearly satisfy 
^ipngahil) = as well as (|18p . Such objects form a 10-dimensional group that may 
be thought of as a generalization of the Poincare group. They have been suggested 
before as useful generators for the linear and angular momenta of extended matter 
distributions [551 UHl HH]- Fixing a hypersurface S that passes through 7 and the 
worldtube of some well-behaved spatially-compact stress-energy distribution T"'', let 

Pa(7,S).l^ + ^5'3b(7,S)B[=bl = / T\i^'^dSa. (24) 

Varying the 10 free parameters here determines the four linear momenta p^ and six 
angular momenta S^^ = S^^^^. Explicit formulae are easily found using ([17]). They 
coincide with standard definitions in flat spacetime (where all t/'k Killing). 
Theorem [2] is easily expanded for these vector fields: 

Corollary 1 Given any Jacobi field f/;^ satisfying C^^gst, — 0, 

C^p^a — C^^Gs — C^^Ga — C^s^CF^a ~ Cr^^H"^^ = 0. (25) 

Again, one argument in each of these equations is assumed to be j. 
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This follows from the well-known identity [531 HI] 

aVa = a^CTa = 2a, (26) 
and its first derivative 

I'a'j'' = (27) 

By definition, a = g3ii{'^)X^X^ /2 is one-half of the squared geodesic distance between 
7 and x. Killing-type Jacobi fields based at 7 therefore drag both arguments of a{x, 7) 
in such a way that distances are preserved. 

It is also possible to identify Jacobi fields that act like homotheties near 7. These 
are distinguished by letting 

5(ab) = ^^^H5ab(7) = cgab(7)- (28) 

As in ([2]), c is an arbitrary constant. For simplicity, the purely Killing components of 
some prospective V'h can be removed by setting A3 — i?[ab] — and c 7^ 0. Substitution 
into P7)) then shows that 

= -ci?"acr' ccr'^. (29) 

This second equality follows from contracting 6^ — — i?"aC^6 with and using (|27p . 

The simplicity of (^5)) is interesting, although perhaps not surprising. It is 
consistent with the interpretation of — cr^ as a "separation vector" between x and 7. 
As has been noted before, a^b ^ in a normal coordinate system. The components 
of cr" would therefore be equal to X". The normal coordinate functions themselves 
act as components of an approximately homothetic vector field. This is unique up 
to a constant factor and the addition of Killing-type Jacobi fields. It generalizes the 
dilations of flat spacetime. 

Generalized Killing tensors of various types can also be generated from Jacobi 
fields. In analogy to (jH) , let 

ICab — ^ipgab- (30) 

These objects exactly satisfy ([S]) at 7, and presumably approximate it near this point. 
It is straightforward to write down other objects which also have this property. For 
example, two (possibly identical) Killing-type Jacobi fields ■0" and ^/j" can be used to 
define 

Kb=^{a^b)- (31) 

This expression clearly generalizes to approximate Killing tensors of any rank. Exact 
second-rank Killing tensors probably exist that cannot be written in either of these 
forms, so it is unclear how useful they are. 

Very near 7, it is possible to approximate the Jacobi fields explicitly. This will be 
especially useful in Sect. [3]below, where a notion of gravitational current is introduced 
with respect to a given vector field. Consider a Taylor expansion of C^gab in powers of 
X^. The first two terms in this series are trivially obtained from (|18p and Better 
approximations involve third and higher derivatives of C^gab hi the limit x —^ "f. The 
lowest order interesting terms can be found from (|A.19p and (jA.24p . Making use of 
(jA.18[) . the final results are that 

^ipgab — a<y^b\j^i)g3h — -^X"^ [C-tP R^chd + 2 ^''^^^f ^acbd )] + 0{X^), (32) 
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and 

2 3 

'^c^iidab — ~-^Cr^a<^ bCr'^cX [i?dc(a -^^i/iSbjf + 5f(a/3^i?b)dc + ^"'^ 

X (25h(a'^^i/>V^i?b)dfc - 5ha.C^Vdi?f(bc)^ " Shb-C^ Vdi?f (ac) )] + O(X^). (33) 

Lie derivatives here are evaluated at 7, so they only involve A^, B^ti, gah, Rated, 
and its first two derivatives. The factors of a^a in front of these equations are 
used as a convenient means for converting tensors at 7 into tensors at x. It is 
perhaps more typical to use parallel propagators g^a for this purpose [24], although 
the aforementioned simplicity of a^a in normal coordinates makes it an attractive 
alternative. There is very little difference at low orders regardless, a^a can be freely 
interchanged with —g^a in ((33|l . This is also possible in (|32p when i?(ab) = 0. 

Approximations like these are not useful over regions where the curvature changes 
significantly, or on length scales approaching the curvature radius. An alternative 
approach is to expand the various bitensors built from a using its definition as an 
integral along a geodesic. Simplifications can often be introduced by ignoring all 
terms nonlinear in the Riemann tensor. A general method for this type of weak-field 
procedure may be found in |23U30| . Specific details involved with expanding the Jacobi 
fields in this way will not be given here. 

3. Symmetries near a worldline 

The Jacobi fields just discussed generalize the idea of a Killing field or more general 
affine collineation in a normal neighborhood of a given point. This is useful for some 
purposes, although it does not have a very clear physical interpretation. The choice 
of origin should presumably correspond to a preferred point, although there are few 
of these that might arise in practice. It is often more useful to base the idea of an 
approximate symmetry off of a given timelike worldline rather than a single point. This 
could correspond to the path of some observer. In some cases, the physical system picks 
out preferred reference frames. A binary star system experiencing no mass transfer 
can admit three center-of-mass frames (rigorously defined in [25[I31[[35] ). for example. 
Two of these are associated with the individual stars, while the third describes the 
system as a whole. There are also preferred observers in most cosmological models. 
Expressing a system's dynamics in terms of quantities associated with these frames has 
an obvious physical interpretation. The distinction between approximate symmetries 
defined with respect to a point versus a worldline is closely analogous to the one 
between Riemann and Fermi normal coordinate systems. 

The concept of an observer here will be taken to mean a timelike worldline T 
together with a set of hypersurfaces E(s) that foliate a surrounding worldtube W. It 
will be assumed that each of these hypersurfaces is a normal neighborhood of the point 
7(5) where it intersects the central worldline. Each of them is therefore formed by a 
collection of radially-emanating geodesies of (usually) finite length. The most typical 
examples would be the past-directed null geodesies or the spacelike set orthogonal to 
7^ = d7^/ds at 7(5). Other choices are possible, however. Regardless, a worldline and 
foliation together will be referred to as an observer's reference frame. 
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3.1. A family of Jacobi fields 

Symmetries adapted to a particular frame can be constructed using a one-parameter 
family of Jacobi fields ip°'{x,"f{s)). Any such family is fixed by specifying Aa(s) and 
Bsbis) as defined in (|19p . An optimal way of connecting initial data between different 
points on T therefore must be found. Before considering this problem, it is first useful 
to collapse the family of Jacobi fields into an ordinary vector field (.""{x). Let t(x) be 
defined so as to identify which leaf of the foliation includes an arbitrary point x in the 
worldtube W. More concisely, it always satisfies x G 'S{t{x)). The assumption that 
each hypersurface is a normal neighborhood of an appropriate point on F implies that 
T is always single- valued. Now set 

eix)^rix,jiT{x))). (34) 

The generalized affine collineations (GACs) to be defined below will be of this form 
for a particular class of famihes ip'^{x,j{s)). 

One potential application for a generalized symmetry constructed using a 
particular frame is in the definition of quantities that might be approximately 
conserved as one moves along F. As an example, consider integrals of conserved 
stress-energy tensors similar to (|24p . One might define the component of momentum 
generated by a of the form ((34| to be 

V^s) = / T\edSa. (35) 

The evolution of this quantity crucially depends on how the parameters ^a(s) and 
-Bab(s) in (fT9ll arc connected along F. It is well-known that is conserved if is 
Killing and no matter flows across the boundary of the worldtube. In this case, initial 
data for the one-parameter family of Jacobi fields must satisfy the Killing transport 
(KT) equations on F: 

DAa/ds =7^Bba (36a) 

DSab/ds = -i?abc^7'Ad. (36&) 

If there exists an exact Killing vector such that A^ = and B^^ = W^Y^ at a 
given s — sq, relations like these will hold for all s. Furthermore, momenta and S^^ 
identified using ([24]) would satisfy 

= (Pa - ^S^'-Rtcd.ilY^ + ^(^ab - 2p[a7b])V = yj^. (37) 

If there were a full complement of ten Killing vectors, all possible versions of this 
expression would together be equivalent to the Papapetrou equations. More generally, 
Papapetrou's result is only an approximation. Any deviations can be understood using 
a general 10-parameter family of possibly approximate isometrics. One might expect 
these corrections to be minimized if A3 and _Bab always satisfy the KT equations even 
when no exact Killing fields exist. 

Definition 1 Let a generalized affine collineation (GAC) SJ^ix) associated with a 
reference frame {F, S} he derived from a family of Jacobi fields via {3^^ . Individual 
elements of the family and their first derivatives satisfy the Killing transport equations 
m^} and (MW on F. 

Although this definition was motivated by the properties of conserved momenta 
in very particular spacetimes, it also arises from much more general (if less physical) 
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arguments. Consider all possible initial data for vectors built from Jacobi fields using 
(|34p . It is reasonable to suppose that any GAC should be exactly affine on F; i.e. 

Va/:5.gbc|r = 0. (38) 

It can also be expected that A3 and Bab fix and its first derivatives on F just as 
they do for -0°. Generalizing (fTOll . let 

A^s) = Chis)), B^'is) = V^eilis)). (39) 

We start with the second of these constraints. Directly differentiating implies 

that 

^^gab = C^dab + 2'0(QVf,)T, (40) 

where the Lie derivative with respect to i''^{x, t) on the right-hand side is understood 
(as usual) to involve only the first argument of this vector field. It will be assumed 
that the foliation is always sufficiently smooth that derivatives of r remain well-defined 
throughout W , and on F in particular. Evaluating ((40|) on the central worldline 
requires knowledge of ■i/'°(7,7)- Coincidence limits like these are commonly denoted 
with brackets. For example, 

[r]{i)^\\m^r{x,-i{s)). (41) 

x^l OS 

The convention of using different fonts for indices referring to x and 7 cannot be 
consistently applied in expressions like this. No confusion should arise, however. 
Limits like (|4ip are easily computed using Synge's rule [531 |H]. In this case, 

= DA7ds-7'^Bb°. (42) 

It follows that ([55)1 holds for all initial data iff (|36ap is satisfied. 

The other KT equation arises from enforcing ((38| . Noting p8|l and (|36aP , it must 
be true that 

[Va^'] = 0. (43) 

[^"] also has to vanish, although this term is equal to — '/''[Vf,?/'"]. Requiring (|43)) is 
therefore sufficient. Using the same type of procedure as in ([42)) shows that 

[VaM = ^Bab/As + Rabc'^YAd. (44) 

Deriving this is straightforward other than noting that ([22| - although mentioned for 
exact affine collineations - also holds for any Jacobi field at its origin. Regardless, the 
conclusion is that the second Killing transport equation ()36b[) ensures that Va£^gbc 
vanishes everywhere on F. 

Noting that the KT equations have the same significance for general affine 
collineations as they do for ordinary Killing fields [15], it easily follows that 

Theorem 3 The class of all generalized affine collineations associated with a given 
reference frame forms a 20- dimensional group in four spacetime dimensions. Every 
GAC satisfies — on T , and all exact affine collineations are members of this 

class. 
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This is closely related to theorem [T] It strongly supports definition [T] and the intuitive 
identification of GACs with approximate symmetries inside W. 

At least in principle, finding GACs associated with a particular reference frame 
is straightforward. Suppose that ^aCso) and i?ab(so) are given as initial data at some 
7o = 7(so)- The goal is then to determine the ^"(x) satisfying l|39)) at the appropriate 
point. This is done by first applying the KT equations to the given parameters along T 
from 70 to ^{t{x)). The geodesic deviation equation ((2T|l is then integrated between 
this latter point and x using the initial conditions (fT9|) . Both of these operations 
simply require finding the solutions to well-behaved ordinary differential equations. 
Alternatively, f could also be obtained using the explicit expression (fTT]) together 
with the KT equations and ([34]). 

Our prescription for generalizing arbitrary affine coUineations may appear 
somewhat awkward. Killing transport equations are being applied along F, while the 
Jacobi equation is used on geodesies intersecting that worldline. These two procedures 
are not as different as they might appear. Trying to use Killing transport everywhere 
would generically lead to inconsistencies. Derivatives of the field expected from the 
KT equations would not usually match the derivatives computed from itself. Only 
the tangential components of these derivatives can be consistently fixed by integrating 
ordinary differential equations along a collection of radial geodesies. Weakening the 
Killing transport equations to take this into account exactly reproduces the geodesic 
deviation equation. This may be seen by rewriting (|2ip as a pair of first order 
differential equations on geodesies connecting x to ^{t{x)). Denote the unit tangent 
vector to one such geodesic by u°-{l). Also set A°- = ijj"- and Bah = a4'b everywhere. 
It is then straightforward to show that 

DiVdZ = u^Bba (45a) 
u"DB,fc/d/ = -RabMu'-Ad. (456) 

The first of these equations has exactly the same form as p6a|) . while the second is 
essentially (j366|) contracted with Killing and Jacobi transport are therefore very 
closely related operations. The latter does not uniquely propagate Bab from given 
initial data, so it is weaker. These remarks also clarify in what sense Jacobi fields or 
GACs approximate Q or ((22)) . 

3.2. Special cases and properties of GACs 

From a physical perspective, momenta like (j35p should be definable even in the absence 
of any exact isometries. This is most conveniently done with a particular class of 
GACs that generalize only the Killing fields. In analogy to the Killing-type Jacobi 
fields discussed in Sect. 12. 2[ suppose that i?(ab) vanishes on at least one point of 
F. It immediately follows from (j366p that it must actually vanish everywhere. The 
Killing- type GACs therefore form a 10-dimensional group of vector fields satisfying 

^feffablr = 0, (46) 

as well as ([55)) . They may be thought of as generalizing the Poincare symmetries of 
flat spacetime near a given observer. 

It also possible to single out GACs that are approximately homothctic in the 
sense that 



(47) 



Approximate spacetime symmetries 



13 



This requires setting -B(ab) = C5ab- While always possible, the remaining components 
of the initial data cannot be explicitly solved for except in the case when F is a 
geodesic. It then self-consistent to choose i3[ab] = 0. The obvious way of doing this is 
to normalize 7^ to unity and set 

A, = c{s - s)% (48) 

for some constant s. It is easily verified that the given parameters satisfy the KT 
equations. One homothetic-type GAG associated with an afhnely parameterized 
geodesic therefore has the form 

eS=c[a'^ + (r-s)i/^a\7^]. (49) 

As usual, Killing- type Jacobi fields may be added to this without spoiling ((47|) . It 
should also be emphasized that homothetic-type GAGs are not restricted to geodesic 
frames. This is just the case where closed-form solutions of the KT equations can be 
obtained by inspection. 

Many of the properties derived for Jacobi fields in Sect. [2] and the appendix can 
be carried over at least partially for the GAGs. For example, 

£5(7= = C^a^ = (50) 

if the arguments are of the form {x, 7(r(x))) and and -0" are related via This 
may be interpreted as stating that spatial Fermi coordinates are preserved under flows 
generated by Since the hypersurfaces E(s) can be described as a set of geodesies 
intersecting j{t{x)), it will always be true that a'^'-V aT — 0. This means that 

Theorem 4 Given a general GAG C^a^ = C^a"^ = when the arguments of these 
equations are as in 150)) . Killing-type GAGs also satisfy C^^^a = C^^^a^ = with 
the same restriction. 

The identity (|X7)) serves to constrain Lie derivatives of the metric with respect 
to Jacobi fields. A direct analog of this equation for an arbitrary GAG would involve 
an additional term. Despite this, contracting the result with cr^ leads to the simple 
conclusion 

aV^/:^5,fc = 2aV'^B(ab). (51) 

"Purely radial" components of C^gab therefore vanish for Killing-type GAGs. 

Many other results can be carried over in similar ways. One that is of particular 
interest is the behavior of C^gab or WaC^gbc near F. As with the Jacobi fields, it is 
possible to see how close a GAG comes to being affine as its reference worldline is 
approached. Analogs of ([5^ and ([55]) may be obtained using expansions like and 
the identity (jA.12|l . Simplifying terms with the Killing transport equations, the lowest 
order correction to ([32)1 is 

~ ^a^a%V(ar5b)f7'^'^'A^hcd' + 0(^3). (52) 
Similarly, the first interesting change to ([33|) has the form 

Va£({-^)5f)c ^ -^Cr^aCT%Cr^c^''y [5h(bVc)T(53 - 7' VaT)£^i?f (Id) 

+ ^VaT(i?df(b'A.9c)h - <?h(bAi?c)fd')] + 0{X^). (53) 

As before, the magnitudes of these terms depend on how close ^° is to being a 
symmetry of the Riemann tensor on the observer's worldline. 
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4. Mechanics and conservation laws 

It has already been remarked that one of the main apphcations of exact symmetries in 
physics is to the formulation of conservation laws. These take several forms. Perhaps 
the most basic are those associated with a spacetime's geodesies. It is sometimes 
possible for such curves to be at least partially parameterized by a number of constants 
associated with geometric symmetries. More interestingly, conservation laws can also 
be associated with extended matter distributions. Assuming only that stress-energy 
tensors satisfy 

VaT"'' = (54) 

tends to lead to the definition of slowly-varying parameters like those discussed in 
connection with (j24p and (|35p . The situation becomes much more interesting in full 
general relativity. Einstein's equation implies stress-energy conservation, although it 
also connects symmetries of the geometry to those of the matter distribution (and vice 
versa). This allows the introduction of exact conservation laws in arbitrary spacetimes. 

4.I. Geodesies 

It is well-known that any Killing fields that may exist provide first integrals of the 
geodesic equation. These can be used both to derive and parameterize the geodesies 
of a given spacetime. While less commonly discussed, similar quantities can also 
be associated with other kinds of symmetries \W[ [33] . Unlike in the Killing vector 
case, the presence of more general coUineations sometimes implies the existence of 
interesting conserved quantities that are not linear in the geodesic's four-velocity. 
The curve's affine parameter can also appear explicitly. As a direct calculation will 
easily verify, two constants associated with an exact affine coUineation are 

Ci = ry'CY^gab (55a) 

C2^yaYl-^lCi. {55b) 

These quantities remain fixed along any affinely-parameterized geodesic y(l). The 
first becomes degenerate if is Killing. C2 then reduces to the standard conserved 
quantity associated with a Killing field. Other constants can sometimes be written 
down by combining Y^ with an exact Killing tensor [111 132] ■ Such constructions will 
not be discussed here. 

Consider instead expressions like those just given with Y^{x) replaced by some 
Jacobi field ip°'{x,j). We then have 

Ci = ^ = -y^ = ry'rViaC^g^cy (56) 

The tangent vectors are proportional to a°'{y,'^) for the special case of geodesies 
passing through 7. It then follows from ()A.10p that both Ci and C2 remain conserved 
along all such trajectories. Each Jacobi field generates exact geodesic constants in this 
way. In terms of the initial data A3 and _Bab, these have the values 

Ci =22/^y^J3(3b), C2=yM3 (57) 

when the parameter I is chosen to vanish at 7. It is clear from this that ^[ab] is 
irrelevant. Multiple Jacobi fields may therefore generate the same constants on a 
particular curve. 
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Exact affine collineations generalize these results by also applying to non-radial 
geodesies. Expansions like ([55)1 can be used to derive how close general Jacobi fields 
come to this ideal. To lowest nonvanishing order, 

Ci ~ -^(2/V^)(j)Vfc)(yV=c)^'i?db//:^.9af + 0(^2). (58) 

This term vanishes if -B(3ti) = 0, so the Killing-type Jacobi fields usually provide more 
accurate "conservation laws" for arbitrary geodesies near 7. In these cases, Ci (1) scales 
like {X/TZ)"^ /TZ, where 7?, is a curvature radius. This is a worst-case estimate. Ci and 
C2 will probably vary much more slowly if there is a physical sense in which the system 
is approximately symmetric. Rates of change for Ci or C2 can also be constrained 
using identities like (jA.Sp . This effectively restricts how much these parameters can 
vary as a geodesic moves away from 7. More of their changes tend to occur as y{l) 
moves across rather than with the radial geodesies. 

Parameters like Ci and C2 can also be defined with respect to a GAC These 
should remain approximately conserved for geodesies near an observer's worldline 
rather than geodesies near a point. They are exact constants for curves passing 
through r along the reference foliation. This can be seen by using (|40)) to show 
that 

ry'r^^aC^gtc = ry'r^aC^gbc (59) 

when y°'V aT = 0. Values of Ci and C2 here are the same as in ([57)1 if the quantities 
in that equation are evaluated at the point where y intersects F. This might be useful 
in a coordinate system constructed from some collection of GACs. Alternatively, it 
can be viewed as a generalization of standard results near a given observer. 

Further methods of parameterizing geodesies may be found by generalizing the 
constants associated with higher-rank Killing tensors. Given any exact Killing tensor 
Kai-ar, = -f'^(Qi---a„), t^c scalar 

CK = K,^-aJ'''---r'' (60) 

is conserved along all geodesies y(l)- An analog of this quantity for the approximate 
Killing tensor ([5U|) is exactly Ci defined above. Something more interesting can be 
generated by substituting (PT]) into (pD|) . In the second-rank case, one may define 

Ck = (j)aV'K)(ybViK) (61) 

for some Killing- type Jacobi fields ■(/'k ^^"^ V'k- This is easily generalized to involve 
an arbitrary number of products, although it is only interesting to consider linearly 
independent collections of Jacobi fields. The maximum number of useful products is 
therefore ten. All of these can be generated just from individual terms of the form 
i'aV'K- These are interpreted as approximate constants associated with objects that 
are nearly first-rank Killing tensors (i.e. Killing vectors). Everything of interest here 
can therefore be derived from the behavior of 

C3=2/aV'K=^^2 + ^^Ci. (62) 

Although this depends only on the two approximate constants defined before, it may 
be interpreted as an additional useful parameter. It has the interesting property that 

C-s = \ci. (63) 

Time derivatives do not appear on the right hand side. Consider the special case of 
a geodesic that passes through 7. Since the Jacobi field was assumed to be Killing at 
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its origin, Ci = 0. This is true everywhere, so C3 also reraains fixed along the entire 
geodesic. It actually coincides with C2 in this case. 

Differences arise when considering non-radial geodesies. It was remarked above 
that there was a sense in which Ci and C2 only varied due to non-radial components 
of y". This type of statement can be made much more precise for C3. Let 

r{l)^u\\{l)a-iyil),j)+ulil), (64) 
where uj[^cr''l = 0. Now (|A.7j) . (|55a|) . and §3^ show that 

C's = ^u^u^C^gab- (65) 

This result is exact for all geodesies y{l). There are many cases where u° becomes 
vanishingly small as I — s- ±00 (when the geodesic exists for these parameter values), 
so (|65p provides a strong restriction on how much C3 can vary in any given situation. 
Very near 7, can be used to show that 

C3 ~ -i«a^)(4a^)X^X^£^i^3cbd + 0(^3). (66) 

The lowest order contributions here scale like {X /TVf' /TZ. This is similar to the result 
expected for Ci when computed using a Killing-type Jacobi field. 



^.2. Extended matter distributions 

From a physical perspective, it is often more interesting to consider possibly 
approximate integrals of the equations of motion describing an extended matter 
distribution rather than a pointlike test particle. Suppose that this matter is modeled 
by a conserved stress-energy tensor T"**. Contracting it with any exact Killing vectors 
that may exist yields conserved currents. These are equivalent to some subset of the 
typical laws of linear and angular momentum conservation known in flat spacetime. 
More generally, ([M)) shows that 

v,(r%y^) - h^'^^'Cygab- (67) 

This holds for any vector field y, although it is convenient to assume that it is a 
Killing- type GAC. The source term on the right-hand side may then be considered 
small near T . This therefore serves as an approximate conservation law. As long as 
there is no matter flow through 9S, quantities like (j35|) might be expected to vary 
slowly in time. 

Stress-energy conservation can be seen as a consequence of the diffeomorphism 
invariance of a system's underlying action. Constructions that are based on it are 
therefore useful in many theories of gravity besides general relativity. They can also 
hold for test bodies in fixed background geometries. This generality is quite restrictive. 
Much more can be said if the full Einstein equation is assumed to hold. Symmetries 
in the geometry are then related to symmetries in the matter fields. The presence of 
an exact Killing field that also satisfies CynTab = allows many interesting results 
to be proven regarding the momenta Pa and S^b defined in (|24p . For example, the net 
force and torque on a body can be written explicitly in terms of the Killing field and 
its first derivative at a point. If is timelike, a body's center-of-mass can be shown 
to follow one of its orbits. Furthermore, mass centers always lie on the central axis of 
cylindrically symmetric spacetimes [12[ 134) . 
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Momenta defined in terms of T°'b^^ are useful for many purposes, altliougli they 
are not conserved in the absence of exact KiUing fields. Determining how they vary 
over time requires detailed knowledge of a body's internal structure. Alternative 
definitions for the linear and angular momenta of an extended body arise when using 
the full Einstein equation rather than just ((54|) . Taking the trace of Ricci's identity 
and rearranging terms shows that 

R^y" = - 9''''g''')^bCYgcd + VfeVi'^y"!. (68) 

This holds for any vector field Note that the second term on the right-hand side 
is always conserved. It follows that 

V„[2(T%r'' - hY-T\)+3^] = 0, (69) 
where the "gravitational current" jy associated with y° has been defined by 

3y = ^(ff^'ff'" - 5""5'")Vb/:y5cd. (70) 

It clearly vanishes if is an exact affine coUineation. This is not the only case where 
the current's contribution to (|69p disappears. Using the Bianchi identity, 

VajJ = -^g'^^'CvRab- (71) 

Any vector field satisfying g'^^'CyRab = will therefore generate conserved matter 
currents involving only T'^^Y'' — FT^{,/2. That such "contracted Ricci collineations" 
generate conservation laws for matter distributions has been noted before in [5D1 [H] . 

The viewpoint here will be to apply ((68|) and ((69)) with Y°- replaced by some 
approximate symmetry. 

Definition 2 Fix some family of Jacobi fields 'ip°'{x,"/{s)) that generates a GAG via 
{34^ . Define the generalized Komar momentum associated with these fields by 

n(^) = Vl^'ld^afc. (72) 

As the name suggests, this has the same form and interpretation as a Komar integral. 
It is convenient to assume that s is a fixed parameter for the purpose of evaluating the 
derivative in (|7^ . Directly using a GAC in place of -0" would add a dependence on 
the reference frame. Note that no such distinctions had to be made for the defined 
in ([55)1 . Applying Stokes' theorem together with ([55)1 and ([70)1 shows that 

= ^ [2(T%V' - \rT\) + ]'i\dSa. (73) 

There is a well-defined sense in which changes in this quantity are determined by a 
combination of "gravitational wave" and matter fluxes across the boundary dY,. In 
this interpretation, the amount of carried away from a system via gravitational 
waves vanishes if the GAC associated with ij}'^ is an afhne collineation. 

The 20-parameter family of scalars is intended to define the linear and angular 
momenta of an extended body. This is at least the interpretation for the 10-parameter 
subset satisfying i?(ab) — 0. As in it is possible to write these momenta in the 
more conventional form of tensor fields on T. Let 

V*^^p:A^ + \st^B^'' (74) 



Approximate spacetime symmetries 



18 



for all and i?ab- As written, the generalized angular momentum 5**1^ needn't have 
any particular index symmetries. Non-Killing Jacobi fields generate the symmetric 
components of this tensor, although such generality isn't necessary. Varying among 
all combinations of initial data completely recovers p* and S*^^. Direct expressions can 
also be obtained with the use of ([17]). Continuing to work with the less explicit form 
(|74p. rates of change of the tensor momenta may easily be extracted from V^. Using 
the KT equations (|36a|) and (j366[l . 

n = {?: ~ \s^.R'',.i')A^ + \{s:, + 2%pI)b^\ (75) 

This is closely analogous to (1571) . The left-hand side is parameterized entirely by A^ 
and Bab, so varying these quantities determines all of the corrections to the Papapetrou 
equations. 

Two definitions have now been suggested for the momenta of an extended body. 
The first - summarized by ([24]) and ([35]) - is closely related to the one given by Dixon 
PSI UHl ISH] • It is well-adapted to the construction of multipole moments for T°-^ that 
intrinsically take into account stress-energy conservation. Mass centers defined from 
these momenta are known to have most of the properties one might expect [I2ll35] . The 
boundary of the worldtube W isn't important as long as it lies outside of the matter 
distribution under discussion. There is no vacuum momentum under this definition. 
Unfortunately, there does not appear to be any exact analog of Gauss' law either. The 
momenta of a matter distribution (and changes to it) must be computed by integrating 
over 3-volumes. The generalized Komar integrals defined by (j72[) have complementary 
characteristics. Their main advantage is in having a direct interpretation analogous 
to Gauss' law. Changes in the component of an isolated body's momentum generated 
by only depends on the gravitational flux passing through the surface dW . The 
mass and angular momenta expected from this definition also agree with commonly- 
accepted notions at least in appropriately symmetric spacetimes. It is potentially 
problematic that includes what is effectively a vacuum energy. These scalars 
usually depend on the spatial extent of W even when its boundary lies far outside of 
any matter distribution. This can make it difficult to neatly separate the properties 
of disjoint matter distributions, although similar situations are found even in ordinary 
electromagnetism. It might be conceptually simpler to extend dW to infinity, although 
it is unlikely that all of the bitensors used here would remain well-defined. There might 
also be convergence problems. Related concepts presented in |36j could be more useful 
for defining momenta over very large regions. 

Some insight into the behavior of the defined here can be gained by computing 
it for very small spheres. To be specific, let C{r,s) be the closed 2-surface on S(s) 
satisfying X^X^ ~ r^ for some r > 0. This is effectively a sphere of proper radius r 
centered at the point 7(5). It follows from ([38]) and (ffO]) that is negligible for small 
radii in the presence of matter. The momentum inside C is approximately 

F| ~ ^r'{T\A'' - Ia't\)\/,t + 0{r*). (76) 

It is perhaps more interesting to consider regions that are locally devoid of matter. 
These can be understood from the behavior of the gravitational current. Its 
approximate behavior near F is easily calculated from (|55|) and ([70]) . To lowest 
nontrivial order, 

4 c -l.a\X^[C^,R\ + ^{frR\ + 2i?b^=')/:v..gcd] + OiX'). (77) 
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This vanishes in vacuum for Kilhng-type Jacobi fields. A httle more calculation finds 
the same conclusion at order as well. This contrasts sharply with other quasilocal 
notions of vacuum momentum in general relativity. An extensive review of these 
concepts may be found in [37 . As remarked there, the energy contained in small 
spheres has been calculated using several different definitions. The generic result is 
that it is proportional to the Bel-Robinson tensor, and scales like r^. If this were 
true for the definition suggested here, terms quadratic in the curvature would appear 
at order in the current. These are not found. Vacuum momenta should really 
only be associated with Killing- type symmetries, so the relevant currents defined here 
decrease at least as fast as as r — » 0. Other definitions in the literature find more 
energy in very small spheres. It is not clear how to interpret this, although it might 
have interesting consequences for the use of near zones and related concepts connected 
to the mechanics of compact bodies. 



5. An example: gravitational plane waves 

The discussion so far has mainly focused on the behavior of generalized symmetries 
near the point or worldline used to construct them. With the exception of general 
identities like (jA.7|l , it has not been clear what happens to these vector fields far away 
from their origins. It is therefore useful to consider an example. Given pT|) . the 
Jacobi fields can all be calculated simply by differentiating the world function. This 
makes it convenient to consider spacetimes where a is known exactly. Essentially the 
only examples of this type are pp-waves or assorted cosmological models (see [35] and 
references cited therein). 

In the interest of understanding the generalized Komar momenta (|72p in a vacuum 
spacetime, only pp-waves will be considered here. Coordinates may be introduced such 
that the metric satisfies 

ds^ = -2dudv + a{u)dx^ + b{u)dy^ . (78) 

It can then be shown that one-half of the geodesic distance between points with 
coordinates (u, v, x, y) and (u, v, x, y) is given by [39l |40] 

a = i[a(u, u){x - xf + /3(u, u){y ~ yf] ~{u~ u){v ~ v), (79) 

where 

u — u u — u 

a '^{w)dw b ^{w)dw 

It is clear by inspection that d/dx, d/dy, and d/dv are exact Killing vectors. They 
are not the only ones. All nontrivial spacetimes in this class admit between five and 
seven linearly-independent Killing fields. 

Most pp-waves are effectively "null dust" solutions of Einstein's equation, 
although there are vacuum examples as well. One of these is given by 

a{u) = cos^{Xu); b{u) — cosh^(Au), (81) 

where we assume that |Au| < tt/2 in order to avoid the two coordinate singularities. 
This represents a simple plane-fronted gravitational wave with amplitude A. The only 
non- vanishing components of the curvature are 

r' X _ y _ \1 lacy\ 
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It is trivial to modify this spacetime to be flat for (say) u < 41j, although impulsive 
waves of this type will not be discussed here. 

Given ([75)l . it is straightforward to explicitly compute the Jacobi propagators 
H°'ba^s and H'^^ah using (fH|) . This can be done for any choice of a{u) and 
b{u), although we will specialize to the case defined by ([8T|l . The results are not 
particularly enlightening to write down in detail, although they have some interesting 
consequences. First, all Jacobi fields are found to be exactly Killing if their first 
derivatives vanish at the base point (denoted as usual with sans-serif font). The 
Jacobi fields i/^jxCv] and i?"[y(Tv] are also Killing. This identifies six independent 
Killing fields. It also implies that the linear gravitational momentum p* defined in 
(|74p must vanish. This can be taken to imply (unsurprisingly) that the gravitational 
wave has zero rest mass: = 0. 

There remain four non-Killing Jacobi fields with skew-symmetric Bab- The one 
associated with spatial rotations in the x ~ y plane is relatively simple to write down 
when u = 0: 

/ tanh(Au) 



d-x^^' 



x) 



/ tanh(AM) 



\ tan(AM) 
- tan(Au) 



d_ 
dy 
d 
dv 



IT.- (83) 



tan(Au) tanh(Aw) 

This clearly reduces to its expected form when u —^ u. The degree to which it succeeds 
in being a genuine Killing field may be estimated by noting that 



\^i>,.y9ab\ 



cos(2Au) -|- cosh(2AM) — 2 



-(Am)* + 0(u^) 



(84) 



V2 sin(AM) sinh(Ait) J 9" 

when u = 0. The quadratic growth estimate here is typically quite good even near 
the coordinate singularities. There is little qualitative change in the nature of this 
expression if u 7^ 0. More interestingly, the gravitational current (jTO]) associated with 
ip^^ always vanishes. This suggests that an observer would not be compelled to ascribe 
any xy component of angular momentum to gravitational waves with the given form. 
iJ°[ijCrv] has similar properties. It satisfies an equation almost identical to and 
the gravitational current generated by it always vanishes. 

More interesting are the remaining two Killing- type Jacobi fields tp^^ = 2iJ''jxO'u] 
and i/jy^ = 2_ff°[y(Tj]. Specializing again to the case u = 0, 

\C^^^gab\' - \C^,.9ab\' - IX*U'{[{X - X)2 + {y - y)^] + l{y - v)u) + 0{u*). (85) 

Unlike the expansion in (|84[) . this approximation fails long before |Au| Tr/2. In 
general, the two magnitudes on the left-hand side have distinct behaviors that strongly 
depend on u. Oscillations generically arise as u is varied, for example. See Fig. 1. 

There are gravitational currents associated with both i/j^^ and V'yu • Applying ([70)) 
and expanding near u = 0, 
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+ 0{u^). 
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Figure 1. Plots of \C^^^gab\'^ for a gravitational plane wave described by I I78I I 
and (HU. The origin is assumed to be at u = 0. Both solid curves assume that 
— V = 0. The dashed ones use \(v — v) = 1/2 instead. Both thicker curves set 
X — X = and \{y — y) = 1/4. The thinner ones use \{x — x) = 1 and y — y = 0. 
Plots for I'C^y^'JabP look very similar unless u 7^ 0. 



These expansions are qualitatively accurate throughout the region of interest. It is 
now clear from (j73p and (j74p that the only non-vanishing gravitational momenta 
(associated with Killing-type GACs) are 



The rates at which these quantities change depends on the relevant fluxes through 
95]. Regardless, the magnitude of the angular momentum tensor always vanishes. 
Intuitively, these statements might be taken to mean that the gravitational wave has 
a "mass dipole moment" equal to its one non- vanishing component of ordinary angular 
momentum. 

The results here could be straightforwardly extended to much more general pp- 
wave (and other) spacetimes. The most interesting point is perhaps the calculation 
of explicit gravitational currents in a vacuum spacetime. In these cases, the general 
results obtained in Sect. 14.21 only state that will decrease no slower than as 
X ^ 0. The example here scales like . Although this conclusion would probably 
not be preserved in more complicated spacetimes, it shows that momenta not arising 
from stress-energy tensors can sometimes be ignored in remarkably large regions. 

6. Conclusions 

Two different notions of approximate affine collineations have been introduced. 
One has the physical interpretation of capturing symmetry principles in a normal 
neighborhood of a point, while the other is adapted to the measurements of a particular 
observer. Flows generated by both of these objects leave a^{x,^) invariant. This has 
the simple interpretation that Jacobi fields preserve Riemann normal coordinates. 
GACs do the same for the spatial components of Fermi normal coordinate systems. 
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These objects always exist, and each forms a 20-diniensional group. Individual 
elements may be interpreted using the values of the field and its first derivatives 
at the appropriate base point. The only caveat to this is that a GAC which might 
initially appear to be purely translational could slowly acquire some rotational and 
boost-type components. This mixing is essential in order to ensure that the fields 
nearly satisfy ([3|) near the observer's worldline. 

The relevance of these definitions ultimately lies in their applications. The 
approximate symmetries introduced here have been used to write down analogs 
of the typical conservation laws applying to geodesies in spacetimes admitting 
affine collineations. Some of the resulting parameters are exact constants of 
motion, while others are only expected to vary slowly near the preferred point or 
observer. Regardless, they may be used to classify and derive geodesies in certain 
regions. Similar results have also been discussed in connection with extended matter 
distributions. This led to natural notions for the linear and angular momenta of 
a spacetime volume as viewed in a particular frame. There seems to be some 
disagreement with other quasilocal notions of gravitational momenta, so it is not clear 
how the definition here should be interpreted. It is unknown if it has any positivity 
or related characteristics. 

Concepts discussed here might also be applied to simplify perturbation theory 
off of some background geometry possessing an exact affine collineation. It could 
be useful, for example, to uniquely construct GACs with respect to a center-of-mass 
worldline that coincides with exact timelikc or axial Killing fields in the unperturbed 
geometry. Center-of-mass trajectories might also be estimated using notions of 
approximate stationarity. More concretely, an analysis of the quantities defined 
in can be shown to provide significant insights into the effects of self-forces and 
self-torques on isolated bodies. Details are presented elsewhere [42] . 
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Appendix: Properties of Jacobi fields 

The Jacobi fields and generalized affine collineations satisfy a number of useful 
identities that both simplify various calculations and further motivate their 
identification as generators of approximate symmetries. The most basic of these results 
is probably (fT8|) or its analog ([38| . Symmetries defined with respect to a given point are 
exactly affine at that point. Symmetries defined with respect to a given worldline are 
exactly affine on that worldline. Nothing is said about the behavior of these fields away 
from these regions. It would be much more interesting if constraints could be placed 
on V aC^pghc throughout the volume N where it is defined. This is indeed possible. 
Jacobi fields generate diffeomorphisms that exactly preserve certain geometric objects 
in all locally well-behaved spacetimes. This fact can be used to derive some properties 
of C^gab and its derivatives even when a; ^ 7. Somewhat weaker comments can also 
be made regarding each GAC Beyond this, approximations for these vector fields 
can be generically obtained near their origins. A few results of this type have been 
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derived before for Killing-type Jacobi fields [551 121] • These will be generalized here 
while also introducing a number of apparently new results. 

The fundamental starting point is the observation that a^(x,j) must remain 
unchanged if both of its arguments are perturbed with a Jacobi field defined about 7 
[29] , More concisely, 

C^a'ix, 7) = + A''(j\ - a^St^ = 0. (A.l) 

This means that Riemann normal coordinate grids are left invariant under infinitesimal 
transformations generated by ip"' . It can be verified using (fT4|) and p^ . The 
solutions to (|A.ip are seen to all have the form ([TT]) derived for the Jacobi fields. 
The undifferentiated world function is left unchanged only under the action of Killing- 
type Jacobi fields. More generally, 

£^(t(x,7) = fTV'^B(3b). (A.2) 

In the special case of a generalized homothetic vector field satisfying (^5]) . this reduces 
to the particularly simple form 

£^jjf7 = 2c(T. (A. 3) 

Returning to the general case, several additional identities may be derived from 
(jA.ip and (|A.2p by differentiation. It is clear that a single covariant derivative with 
respect to x will commute with the Lie derivatives in both of these equations. Using 
this fact shows that 

C^cj\ = C^,H\ = 0, (A.4) 

and 

C^,(7a = 2a^a^aB(^^b)- (A. 5) 

It then follows from and (|A.4p that 

C^a" = 0. (A.6) 
Comparing this to what would be expected from (jA.5|) produces the very useful identity 

a^C^Qab = 2(TV''f,B(at,). (A. 7) 

This result is central to the interpretation of Jacobi fields as generalized symmetry 
generators. It applies most clearly to Killing- type vectors V'k- The right-hand side 
then vanishes, providing a strong restriction on how much these fields can fail to satisfy 
Killing's equation away from 7. 

More generally, there should exist an analog of (|A.7p that constrains V aC^gbc for 
all Jacobi fields. Expressions of this form are easily obtained by directly differentiating 
(jA.7|) . Applying identities like ([?71) to the result then shows that 

aV''[2V(a/:^fff,)c - ^cL^gab] = (A.8) 

and 

a'^a^VaC^gbc = 2a= - a^a^6)B(3b) • (A.9) 

The first of these equations ensures that that there is always some sense in which 
every Jacobi field is approximately affine. It is actually identical to ([51]) . One trivial 
consequence is that 

aVVVa£v5bc = 0. (A.IO) 
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It means that the "purely radial" component of ^ holds everywhere that -0" is defined. 
Similar statements can also be made using (jA.9p . The right-hand side of that equation 
vanishes for Killing or homothetic-type Jacobi fields, although interesting remarks 
can be made even in the general case. Temporarily consider points that are spacelike- 
separated from 7, so fj > 0. It is then clear that in normal coordinates, the appropriate 
components of VaC^gtc fall off roughly like l/v^ far away from the origin. These 
quantities also have a simple scaling behavior near 7. Given that all third derivatives 
of the world function vanish when x ^ 7, the term 

a\ - cj\a^b = (T'^a'^ab (A.ll) 

in (|A.9|) must go to zero at least as fast as a in this limit. It then follows that the 
portions of Va^ipgbc constrained in that equation should grow no faster than y/a for 
X very near 7. 

A more explicit connection can actually be made between Va^^^jgbc and 
derivatives of the world function. Several remarks have already been made about 
the interpretation of various results here in terms of normal coordinate systems. Since 
([3| is essentially the condition that leaves the Levi-Civita connection invariant, one 
might expect an appropriate result to arise from studying the Christoffel symbols in 
normal coordinates. These are related to the bitensor a^ab- Lie derivatives of this 
quantity can be computed by directly differentiating (jA.4p . Simplifying the result 
requires the relation 

VfcVaV'c = Rcabdi''^ + ^[2V(a/:v5'&)c " c^^Qab]- (A. 12) 

This is obtained by repeatedly applying Ricci's identity, and actually applies for any 
vector field. It is a generalization of and the well-known expression for second 
derivatives of a Killing field. Combining ((M]), (jA.4[) . and (|A.12[) shows that 

V aC^gbc — 2i?(c|a|'C^Cr^6)a- (A. 13) 

This is essentially the expected result. 

There is at least one more quantity derived from the world function that has a 
simple interpretation in normal coordinates. That is the van Vleck determinant 

This two-point scalar arises as part of the volume element in these charts. It also 
appears in connection with the focusing of geodesic congruences and in Green functions 
associated with common wave equations [24]. It is clear from the definition (jA.14[) 
that 

VblnA = -iJ^a^b. (A.15) 

Derivatives with respect to 7 arc identical if the index h in this equation is changed 
to a b. Simplifying with (|A.4p then shows that 

Z:^ In A = - VaV-" = B\ - ^g''''C^gab. (A.16) 

The results obtained so far apply everywhere that the generalized symmetries are 
defined. It is sometimes more useful to derive explicit approximations for the Jacobi 
fields. These can be found in various ways. Perhaps the most obvious procedure is to 
evaluate (|17p in a normal coordinate system with origin 7. The second term involving 
B^^ then reduces to a linear combination of the coordinate functions. This is exact. 



Approximate spacetime symmetries 



25 



The factor of fi^ti appearing in the other term does not simphfy in these coordinates. 
It reduces to the identity at the origin, ahhough its behavior further away must 
be computed in some other way. For very small, known Taylor expansions for 
(7^b could be used [H]. There are also integral approximations neglecting all terms 
nonlinear in the curvature [23l [30]. Either of these methods can be used to find 
explicit expressions for ip"' or in appropriate situations, although it is often more 
interesting to directly compute Lie derivatives of the metric with respect to these 
vector fields. The viewpoint here will be to derive a Taylor series for C^Qab centered 
at 7. This expansion will be in powers of the radial vector —a^. The first two terms are 
clear from and pS]) . Going to higher orders requires knowledge of [VaVh/I^gcd], 
[VoVbVc/C^gd/], and so on. The bracket notation here was defined in (|4T|) . In any 
case, limits like these may be computed from (jA.lSp . Results such as ()A.12|) can 
further be used to derive the full third or fourth covariant derivatives of ip'^ in the 
coincidence limit. 

First consider third derivatives. From (|A.12p . V b^^dcd involves [Vf,-ff°a] and 
[Vd-Cj/iCr^fjc]. All third derivatives of the world function vanish at coincidence, so only 
the second term is important here. Some manipulation shows that 

The quantity being differentiated on the right-hand side is known to have the form 

[SUES! [30] 

2 

[o^bcd] = -^R^(bc)d, (A. 18) 

so 

2 

[Va^bC4,gcd] = g (^ab(/^i>.9d)f + .9f(c'Cv.-Rd)abO- (A.19) 

Lie derivatives here have their usual expansions in terms of Ag, Bab, ^abc'', and 
Vf/Jabc"^- The first term clearly vanishes for Killing- type Jacobi fields. The 
other depends on how close comes to being a curvature collineation (satisfying 
^YcRab/ = 0) at 7. 

It is not too difficult to compute one more derivative of C^gab- This is again 
derived from (jA.lSp . Noting that 

[C^C7\,] = C^a'bc] = 0: (A.20) 

one finds 

[VaVf,Vc/:v5d/] = 2[VaVf,£^c^^(<^].9f)h. (A.21) 
The two covariant derivatives here may be commuted inside the Lie derivative. A 
number of remainder terms arise, although these all vanish on account of (|A.12[) and 
(|A.18|) . This leaves an expression dependent only on [cr^bcd/]. As noted in [23] . 

Wabcdf] = V(fi?d)(ab)c + 2^c^^d(ab)f- (A. 22) 

A straightforward application of Synge's rule together with (jA.lSp shows that 

W^bcdf] = gV" 

Combining all of these results. 



W^bcdf] = g V''i?b(df)c - V(fi?d)(bc)''- (A. 23) 



[VaVbVc£^5d/] — g.9h(d^i/'^'^-^f)(ab)c ^ 3hd^i/)V(ai?b)(cf)'^ 

-5hf>Cv,V(ai?b)(cd)'- (A.24) 
Taylor expansions derived from these limits are given by ([32| . ([33]) . ([SS]) . and ([53]) . 
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